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Abstract 

The power of the disconjugacy properties of second-order differential equations of 
Schrodinger type to check the regularity of rationally-extended quantum potentials 
connected with exceptional orthogonal polynomials is illustrated by re-examining the 
extensions of the isotonic oscillator (or radial oscillator) potential derived in feth-order 
super symmetric quantum mechanics or multistep Darboux-Backlund transformation 
method. The function arising in the potential denominator is proved to be a poly- 
nomial with a nonvanishing constant term, whose value is calculated by induction 
over k. The sign of this term being the same as that of the already known highest- 
degree term, the potential denominator has the same sign at both extremities of the 
definition interval, a property that is shared by the seed eigenfunction used in the po- 
tential construction. By virtue of disconjugacy, such a property implies the nodeless 
character of both the eigenfunction and the resulting potential. 
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I INTRODUCTION 



The introduction of the concept of exceptional orthogonal polynomials (EOP) [TJ [2] and 
the construction of the first Xi families of Laguerre- and Jacobi-type EOP in the context 
of Sturm-Liouville theory [3j H] have aroused a lot of interest both in mathematics and in 
physics. On the one hand, the EOP families are novel systems of complete orthogonal poly- 
nomials with respect to some positive-definite measure, generalizing the classical families 
of Hermite, Laguerre, and Jacobi, but differing from these among others by the presence 
of gaps in the degrees n of the polynomials appearing in the sequences (e.g., n = for the 
X\ families). On the other hand, they were shown [5] to provide a keystone for rationally 
extending some well-known (translationally) shape-invariant quantum potentials J6J [8] 
in such a way that shape invariance is preserved. In this context, the usefulness of an 
approach based on supersymmetric quantum mechanics (SUSYQM) or, equivalently, the 
Darboux transformation was pointed out [HI fTTJj . 

An important step was the description of infinitely many shape-invariant potentials as- 
sociated with X m EOP families, where the number of gaps m, called codimension (and 
corresponding here to n = 0, 1, . . . , m — 1), may be arbitrarily large [HJ [12j 1131 E]. The 
existence of two families of X m -Laguerre and X m - Jacobi EOP was also reported (thereby 
extending an observation made for m = 2 in Ref. [10]) and it was later on explained through 
Darboux-Crum transformation JT5J [THJ [T7] . In addition, these families were shown to be ob- 
tainable through two alternative approaches, namely the Darboux-Backlund transformation 
[T5] and the prepotential method [19J. 

Counterparts of the classical family of Hermite polynomials have also been known in 
physical applications since the early 90s EB E21 E31 1211251 ESI EZ] , but these X m EOP are 
of a slightly different nature in the sense that the m gaps in the degree sequence correspond 
to n = 1, 2, . . . , m, instead of n = 0, 1, . . . , m — 1. Since then, similar X m -Laguerre and 
X m - Jacobi EOP have been constructed and form the so-called third families [101 HH [T9] . 

A significant advance was then the introduction of multi-indexed A miim2i ... jmfe EOP 
families and associated rationally-extended potentials through multistep Darboux algebraic 
transformations [2H] . The same concept was developed through several approaches, such as 

2 



the Crum-Adler mechanism [29], higher-order SUSYQM [301 [31], an d multistep Darboux- 
Backlund transformations [32J. 

In such a context, there remain two open essential questions. The first one is to 
know whether the multi-indexed EOP families exhaust all the possibilities of higher- 
codimensional complete orthogonal polynomial systems or, in other words, whether all the 
higher-codimensional complete orthogonal polynomial systems are generated by the appli- 
cation of successive algebraic Darboux transformations (or one of the equivalent methods 
quoted above). This was recently proved to be true for codimension two [33] . 

The second problem has to do with the existence of a well-behaved measure defining an 
EOP system, which is directly connected with the regularity of the associated rationally- 
extended quantum potentials. In the one-indexed X m case, this problem is usually solved 
by showing that the classical polynomial occurring in the denominator has no zero in the 
domain of the variable, which is an easy task since the distribution of zeros of classical 
polynomials is well known [3H [35]. In the multi-indexed X mi m2j mfc case, however, the 
denominator contains a Wronskian of functions written in terms of classical polynomials, 
so that the problem becomes much more difficult to solve. While still feasible for k = 2, 
the study of zeros remains outside scope for most higher k cases, hence the construction of 
multi-indexed EOP is still rather formal. 

A notable progress towards the solution of the regularity question was recently made [32] 
by making use of the disconjugacy properties of the Schrodinger equation for eigenvalues 
below the ground state [361 EH EH]. As a consequence of such properties, the nodeless 
character of the function present in the denominator can be inferred from the signs it takes 
on both ends of the domain. Nevertheless when applied to the rationally-extended isotonic 
oscillator (or radial oscillator) potentials [32], except in the k = 2 case, the equality of signs 
was obtained without exhibiting explicit formulas for the denominator in the asymptotic 
regions. 

The purpose of the present paper is to give such explicit expressions in the general multi- 
indexed case. We then extend some results previously obtained in [301 [31 ], confirming in 
particular a conjecture about the structure of the polynomial factor in the denominator. 
This allows a more direct proof of the regular character of the considered extensions. 
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The paper is organized as follows. After recalling in Sec. II the disco njugacy properties 
of the Schrodinger equation, we apply it to the construction of rationally-extended radial 
oscillator potentials and corresponding X m -Laguerre EOP in first-order SUSYQM. Section 
III summarizes the corresponding multistep construction in fcth-order SUSYQM, where 
each step may be associated with any one of the first two families of X m -Laguerre EOP. 
An improved proof of the regularity of these extensions is then given in Sec. IV. Finally, 
Section V contains the conclusion. 



II RATIONAL EXTENSIONS OF THE ISOTONIC 
POTENTIAL AND X m -LAGUERRE EOP 



A Disconjugacy 

Let us briefly recall some essential elements about the disconjugacy properties of second- 
order differential equations of Schrodinger type, 



To any solution ip{x) of Eq. (12.1 1) we can associate the corresponding Ricatti-Schrodinger 
function w(x) = —tp'(x)/ip(x), which is a solution of 



Equation (12.11) is said to be disconjugated on I C R {V(x) being supposed continuous on 
/) if every solution of this equation has at most one (necessarily simple) zero on I [36, I3"T] . 
For a closed or open interval /, the disconjugacy of Eq. (12.11) is equivalent to the existence 
of solutions that are everywhere non zero on / [36j EZ] • In the following we will consider 
J=]0,+oo[. 

Let us also quote the following result: 
Theorem |36l [37] // there exists a continuously differentiable solution on I of the Riccati 
inequality 



il>"(x)+ (E -V(x))if>(x) = 0. 



(2.1) 



w'(x) + w 2 (x) = V(x) - E. 



(2.2) 



w'(x) + w 2 {x) < V(x) - E, 



(2.3) 



then Eq. k2. 1\) is disconjugated on I . 
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In the case we consider, the Hamiltonian corresponding to Eq. (12. ip and associated 
to Dirichlet boundary conditions on / has an infinite bound-state spectrum {E u , tp u }ueN- 
Consequently, since ipo is nodeless on J, wo(x) is a continuously different iable solution of 
(12. 2p on this domain and condition (12. 3 p is fulfilled for every value of E not greater than Eq. 
The domain of values ] — oo, Eq] of the spectral parameter E will be called the disconjugacy 
sector of Eq. (12. ip . In this sector, except for E = E , ip(x) does not satisfy the required 
Dirichlet boundary conditions to be a bound-state wavefunction and diverges at one or 
both extremities of I. To ensure that this eigenfunction is nodeless, it is sufficient to verify 
that it has the same sign at both extremities of the definition interval. 

B Isotonic oscillator and X m -Laguerre EOP 

The isotonic oscillator potential is defined on the positive half-line x > as 



Vt(x) = i w V + (2.4) 
4 ar 



The spectrum of the associated Schrodinger equation has an infinite number of bound 
states, whose (unnormalized) wavefunctions are given by 

^\x)o^ Vl (z)L^(z), ueN, (2.5) 

where z = ux 2 /2, a = I + 1/2, and 

m (z) = 2 («+V2)/2 e -*/2 j (2.6) 

Lv(z) being a generalized Laguerre polynomial [31]. The associated energies are 

E® =w(2i/ + a + l). (2.7) 
In first-order SUSYQM, the superpartner Hamiltonians are defined as [7j 

#(+) = A^A = —ft + Vl+Hx) - E, 



dx 1 

= AA^ = --£2 + V(~\x) - E, 



{21 



with 

d 



A=£ + W(x), 



(2.9) 



where the superpotential W(x) satisfies the Riccati equations 

=F W\x) + W 2 (x) = V {± \x) - E. 



(2.10) 



W(x) is obtained from a solution <p{x) of the initial Schrddinger equation 

H M (f)(x) = 0, (2.11) 

via 

W(x) = -<f>'(x)/<f>(x). (2.12) 

The new potential V^(x) = V {+ \x) + 2W'{ x) is regular as soon as the "seed" solution 
4>(x) is nodeless on ]0, +oo[. This can be achieved if we choose the factorization energy in 
the disconjugacy sector of H^ + \ that is, when E is not greater than the ground-state energy 
E { +) of VW(x). This is true in particular if we choose as seed solution the ground-state 
wavefunction i/Jq\ x ) = r]i(z) itself. For E < E^ + \ as discussed above, the analysis of the 
sign of <fi{x) at + and +oo is then sufficient to control the regularity of W(x). 

Suppose that we are in this case. Then we have three possibilities (up to a global change 
of sign), namely 

I : 0(0+) = 0+ 0(+oo) = +oo, 

II : 0(0+) = +oo, <p(+oo) = 0+, (2.13) 
III : 0(0+) = +oo, 0(+oo) = +oo. 

In the first two cases, l/<p(x) diverges at one extremity of the positive half-line and, 
although we have 

it is not a bound-state wavefunction of H^~\ Then and turn out to be isospectral. 
If we are looking for seed functions which are polynomials (up to some gauge factor), we 
obtain them from the eigenstates by applying the following discrete symmetries acting on 
the parameters of the initial potential (12. 4 p [T8] 

Ti : (u,a) ->■ (-u,a), 

T u :(u,a)^(u,-a), (2.15) 
r m : (cj,a) -> (-u, -a). 
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By applying the symmetry Tj, i — I, II, III, to ipv, we obtain a seed eigenfunction (p v , 
which satisfies boundary conditions of type i (see (12.131) ). In particular, for i = I or II, we 
obtain 

f T,:^ { l ) (x)^dUx) = y](z)Lt ) (-z). 

(2.16) 



Fj : ^ } (x) -> 0L(^) = xi(z)L%\-z), 
r 11 :^(x)^cf>Ux)^x} 1 (z)L^\z), 



where 



X }(Z) = 2 («+l/2)/V/ 2 , 
X JI( Z ) = z -(a-l/2)/2 e -V2 ) 



for the corresponding energies 



Irn 



-u(a + 2m+ 1), 
-a; fa — 2m — 1). 



(2.17) 



(2.18) 



If 4>} m is always in the disconjugacy sector, for 0}^ to reach it we need to assume the 
condition a > m. Starting from the isotonic potential = Vi'(x), the SUSYQM 

partnership based on seed eigenfunctions (f>j, m or <f)]} m with respectively V = I — 1 and 
/' = / + ! gives then rational extensions of the form [31] 



V(-\x) = V l (x) + V ltnt (x) + C, 



where 



Vi tTat (x) 



-2u 



.(a) 
9m 



+ 2z 



..(a) 



• a) 



(the dot denoting a derivative with respect to z), with in the first case 



T (a-l)/ 



c 



-a; 



and in the second case 



» 



C = w. 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



Due to disconjugacy properties, the regularity of these extensions on the positive half- 
line is ensured for every m — 1, 2, . . . (and a large enough for type II). 

The bound-state wavefunctions of the superpartner potential V^~\x) are obtained by 
applying the operator A (with W (x) = —^, m {x)/<j)\i m {x), i = I, II) on those of 
(^ +) (a;) oc r] V (z)L i u l '\z) with «' = /' + 1/2). We can write 



z/ = 0,l,2,..., 



(2.23) 



where Dm+v{ z ) is a polynomial of degree n = m + u, 



y£\z) = { L «r niz) (2.24) 

L a,m,n( z ) lf 1 = H > 



satisfying the second-order differential equation 

d 2 ( n (a) \ d n (a) n {a) 

Z-— + Q; + 1 — Z — IZ—r-r — + [Z - OLJ i~t + Z j r 

dz 2 \ / dz v V a) ( a ) 

\ ym I ym ym 



Vmlviz) = -vy { m+v{z)- (2.25) 



The set of polynomials y^m+v{ z ) (f = 0, 1, 2, . . . ), defined on the positive half-line, is 
orthogonal and complete with respect to the positive-definite measure of density 

*re-/(g£\z)) 2 . (2.26) 

III LAGUERRE EOP AND REDUCIBLE fc-TH OR- 
DER SUSYQM 

In /cth-order SUSYQM, the first-order differential operators A and A' are replaced by kth- 
order ones A and A^ [22J [231 Ell HOI |JT] . The correspondence between the initial and final 
Hamiltonians and is determined by the intertwining relations 

AH^ = H^A, A ] H^ = H^A ] . (3.27) 

In the reducible case, the operators A and A^ can be factorized into products of first- 
order differential operators as 

A = A ik) ■ ■ ■ A {1 \ (3.28) 

where 

A® = ■$- + W®(x), W {i \x) = -(log0 (4) (x)V, (3.29) 
dx 



the seed eigenfunctions being given by Crum's formula 

M')=<® M- X^-^l (3-30) 
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with W(0i 



\x) denoting the Wronskian of 0i(x), 02 (x), ■ ■ ■ , 4>i{ x )- Here we have 



introduced lower indices 1, 2, . . . , i to specify which seed functions have been used. The 
partner potentials are then linked by the second Crum's formula 



V® (x) = V® (x) - 2-f^ log W(0i, 2 , . . . , <f> k \x). 

dx z 



(3.31) 

Vi'(x) to get an extended 



As in Sec. II, we start from the isotonic potential V^\x) 
potential of the form 

V {2 \x) = VKx) + V hrs ,{x) + C. (3.32) 

The final potential is of type I q ll ~ 9 if the chosen set of seed eigenfunctions includes 
q < k seed functions of type I and k — q seed functions of type II. The order of the seed 
functions being irrelevant, we can assume that 

ftvmS?) = xU^L^i-z), i = l,...,q, 

. ,. (3.33) 
4™» = X$(z)I&*\z), i = q + l,...,k, 

with I' = I + k — 2q (hence a' = a + k — 2q), < mi < ■ ■ • < m g , and < m q+ i < - • • < mjt, 

it has been shown that the 



(3.34) 



(f)i(x) 



these functions being nodeless for a' > sup (mi). In 

i=q+l,...,k 

Wronskian appearing in Crum's formula (13.311) can be written as 

W(0 1 ,0 2 ,...,0 fc | x) = {u>x)« k -Wz^ k ^(ti) 9 {^gP{z), 



where 



with 



( z ) = z -(k- g )(k- q -i) det fW ) 



(3.35) 



r(«) 



^mj-i+li 2 )' 

{m j + l) i - 1 z h - i L^-if\z) 

(a'+i-l). 



vrij —q 



-z 



(rrij + l) q {mj - a' - i + q + 2)i_q_i 



x z 



k-i t (-a'-t+l)/ 
J m,j+q ' 



1 < % < q+ 1,1 <j < q, 
l<i<q+l,q+l<j<k, 
q + 2 < i < k, 1 < j < q, 

q + 2 < i < k,q + l < j < k, 



(3.36) 



and (a) n = a(a + 1) • • • (a + n — 1) denoting the usual Pochhammer symbol 
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Observe that on taking ( 13. 3 1 p and (I3.34p into account, we can express Vi tiat (x) and C 
of Eq. fl3T32|) as 

2" 



.(a) 

-2u < + 2z 



■■(a) / .a) 



(a) (a) 



C = (k — 2q)u. 



In [31], it has been proved that 



(3.37) 



9u*\ z ) = + lower-order terms, 



(3.38) 



with 



A(m,. 



= Eli m< - q(q - l)/2 - (k — q)(k - q - l)/2 + <?(*; - g), 
Cjr = (-l) <T A(mi, . . . , m,)A(m, + i, . . . , m fc )/ (mj . . . m k \), 

° = J2i= q +i m i + <l( k -<l), 
, 7ty) being a Vandermonde determinant of order I, 



(3.39) 



A(n 1 ,...,n 1 ) = j | (nj-rii). 

l<i<j<l 



(3.40) 



As shown in [32], the Wronskian f !3.34j) has no node on the positive half-line, which 
ensures the regularity of the rationally-extended potential V^ 2 \x) on this domain. The 
bound-state wavefunctions of the latter are then given by 



^\x) = ^-y%{z), 



z/ = 0,l,2,..., 



(3.41) 



where y^. u (z) is a polynomial of degree n = // + v satisfying the second-order differential 
equation 



y%{z) = -vy% v (z). (3.42) 



(«) 



The y^u+v ( u = 0) 1> 2, • • • ) polynomials constitute a complete orthogonal family on the 
positive half-line with respect to the positive-definite measure of density 



(3.43) 
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IV POLYNOMIAL CHARACTER OF AND BE- 
HAVIOUR NEAR THE ORIGIN ^ 

In [31] it has been conjectured that g^\z) is a polynomial, which can be directly checked 
in the k = 2 case, but is not obvious for higher k values whenever k — q>2 (see Eq. (I3.35j) ) . 
In the following we will give a complete proof of this assertion. In [32] the regularity of the 
extended potential, which is equivalent to the absence of node of \ has been established 
without giving the exact highest and lowest degrees coefficients of this polynomial. We 
will also determine the value of the lowest-degree coefficient and then verify explicitly the 
nodeless character of the Crum Wronskian. 

First note that from Eq. (13.341) we can write for k — q > 2 

w(0 1; 2 ,...,0 fc i x) = (ux^y^-^ixiyixj^g^iz), 

W(&, | x) = (coxf- 1 ^- 2 ^ z-^- 1 Kxld q (x}}) k - q - 1 9 { ; +1 \^), 

where I' = I + k — 2q (hence a' = a + k — 2q) and 

» = Eli m i ~ Q(q ~ !)/2 -(k-q)(k-q- l)/2 + q(k - q), 
»' = Eti rrii - q(q - l)/2 - (k - q - l)(k - q - 2)/2 + q(k - q - 1) (4.45) 
— fi — rrik + k — 2q — 1 . 

Inserting these expressions in Eq. (13.301) . we arrive at 



(4.44) 



(2a;) (fc - 1)/2 z- (2a, - 2fc+49+1)/4 e-^ /2 ff^ . (4.46) 



(k) 

On the other hand, we can obtain a different expression for <j)\ £ k by using some properties 
of Wronskians. Indeed, from Sylvester's theorem [13] , we can write 

W(0i, . . . Ak-lAk I x) 



and using 

W(f<t> 



W(W{<f> 1 ,...,<j> k -2,<t>k-i |x),W(0i,..., 




W(0i, . . . ,0fe_2 | X) 






(4.47) 


i,...,f(f>k | ar) = / fc W(0i,...,0 fc | ac), 


(4.48) 
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we have 



x W 



that is 



'l,2,. ..,fc 



w 



x 



dk-i) 



,(*-!) 



'fc-2 1 ar) 






fc-2, 0fc a;) 




4>k-2 X) 


-15 < ) 9 l,...,fc-2,fc 


XJ 


-2,fc-l( X ) 



X 



(4.49) 



(4.50) 



From the Wronskian theorem [44J, it follows that the Wronskian in the numerator of 
this expression satisfies the property 



,(*-!) 



,(fc-i) 



%...,k-2,k-V V y l,...,fc-2.fc 



{Eu-i - E k ) I 0i fc - 1 fc ) -2,fc-i(O0i:X 2 , fc (O^, 



(fc-1) 



(4.51) 



because k-2k-i anc ^ fc-2 k correspond to the energies YTi=i E% and Yli=i Ei + Ek, 
respectively. 

As a consequence of the exponential factor in Eq. (14. 46 p . both 0i~'2_2 fc-i(0 anc ^ 
^1 fc-2fc(0 vanish for £ 00, which ensures the vanishing of their Wronskian in the 
same limit. Hence Equation (14.5ip reduces to 

,(fc-i) ,(fc-i) 



,(fc-i) 



rfc-1 



^fc-2 



'l,...,fc-2,fc-l> 'Tl,...,k-2,k 



X 



r+00 

(E k -i-E k ) / 0£: J i ) _ 2ifc _ 1 (e)0 1 !: J i ) - 2 , fc (e)^, (4.52) 

or, on taking Eq. (14.461) into account and defining ( = w£ 2 /2, 

■ -{2u) k ^ 2 {E k ^-E k ) 



VV \ V 9 ! fc-2, fc-D V 9 ! fc-2,fc 



X 



+00 



Ki) ( . («+i) 

^.fc-2g-a'-2 e -c V \Sl9p,' \ 



<# +2) (c))' 



(4.53) 



where 



(4.54) 



p! = fi — m fc _i + k — 2q — 1, 
/i" = /i — m^-i — m/j + 2(fc — 2g) — 3. 
Inserting Eqs. (I4.46P and (I4.53P in the right-hand side of Eq. (I4.50p . we get the searched 
for second expression of <f>\\ k {x), 



dk) 

h,2,...,k 



. x 



-{2u)^ 2 {E k ^-E k ) 

"+00 



x 



X 



C k-2 q -*'-2 e -t g (« + l) (c) ^ + l) (c) / fe (« + 2) 



2 (2fc-4 g -2a'-3)/4 e -,/2^( ? +l) (z)/ ^+2) (z) 



-1 



(4.55) 
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which can also be rewritten as a recursion relation for the g^\z) functions, 



Eu — E, 



fc- 1 r -(k-2g-a'-l) z (a+2) ( s 



c 



k-2q-a!-2 _f (a+1) 



^ +1) (C)4' +1) (C)/(^ +2) (C))^C (4.56) 



On the right-hand side of this relation, and g^ l> correspond to /c' = k — 1 



>+i) 



while for g^, +2 \ k' = k — 2. The polynomial character of g K ^'{z) can be easily obtained by 
induction. This is already verified for k = 1 and k = 2. Assuming that for k' < k, is 
a polynomial with a non-zero constant term, it turns out that in the neighbourhood of the 
origin, the right-hand side of Eq. (14.561) has a well-defined value given by 



(«)/ 



a' 



-9 



(a+l), 



(0)4? +1) (0)/^ +2) (0)(l + O(z)), 



>+2), 



(4.57) 



fc + 2g + l^ 

where we have used E k = —u(a' — 2m k — 1). 

From Eq. (I4.57P we immediately deduce that g/f > (z) cannot have a pole at the origin 
and is then a polynomial with non-zero constant term. 

To completely describe the behaviour of g^ 1 (z) at the origin, we need to determine the 
exact expression of this constant term. For such a purpose, we will proceed by induction 
again. Consider first the pure case l h (k = q), where a' = a — k. From (13.351) and (13.3 
we have 



Lmi ( z) 



(a'+k-l) ( _ x j 



(a'+fe-l)/_ s 
m k -k+l \ A ) 



On using 



we easily get 



9^(0) 



(a + 1) 



(4.58) 



(4.59) 



(a' + fc) mi _ fc+ i ■ ■ ■ (a' + /c) mfc _ fc+ i 
mi! • • • m k \ 

(a' + 1) • • • (a' + k - 1) 
{a' + 2) • • • (a' + k - l)m x 



x 



{a' + 1) • • • {a' + k - 1) 
(a' + 2) • • • (a' + k - l)m k 



mi (mi — !)••• (mi — k + 2) ... m k (m k — !)••• (m k — k + 2) 



(4.60) 
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that is (see Eq. QS3QD ) 



(«' + 1) 



mi 



{a' + k) 



mi! • • • m k \ 

(a - k + l) mi • • • (a) mk -k+i 



A(mi, . . .,m k ) 



A(mi, . . . ,m k ). 



(4.61) 



mil ■ ■ ■ m k \ 

Consider now the first mixed case (k = q + 1), where a' = a — k + 2 [31]. We 

then get 



(a'+k-2) ( _ x 
(a'+k-l), x 



^^ 1 fc _ fc 2 l 2 (-2) (m k + l) k ^ 2 zL K m * +k 



(-a'-jfc+2). 



(a'+fc-l) 
mfc_!— fc+1 ^ 



-z) (m fe + l) fc _i-L 



(-a'-fc+l), 
m^+fc— 1 \ 



(4.62) 



which leads to 



9i a \0) = (m k + l) k _ 1 L^- k+ 1 1 \0) 



L 



(a'+k-2). 



-(«') 



,(0) 



-/, +2^(0) ••• ^mfe_i-A:+2(0) 
The right-hand member determinant is exactly the constant term of the g polynomial 
in the pure case l k ~ 1 and we deduce from the previous result that 

[a' + l) mi •■■(a' + k- l) mfc _ 1 _ik+2 



(4.63) 



#(0) 



mi! • ■ -m fc _i! 



: A(mi, . . . ,m fc _i) 



or 



,(o 



)(0) 



1) 



x (m fc + l) fc -i(-l)^ +fc - l ( f ™i mk+ *-\ 

(m k + k- 1)! 

-i ( a ' + • • • («' + fc ~ l)m k _ 1 -k+2(a' ~ rn k ) mk+k _i 

mj • ■ ■ m k \ 



(4.64) 



x A (mi, . . . ,m k -i) 

>m fc +fc-l ( Q ~ fc + 3 )"^i •••(« + l)m fc _ 1 -ife+2(Q! ~ fc + 2 - m fc ) mfc+fc _i 



■1)' 



x A(mi, . . . ,m fe _i). 



For the general mixed case PII fc 9 , we now make the following hypothesis 



mi! • • • m k \ 



(4.65) 



si a) (o) 



CT A(mi, . . .,m q )A(m q+1 , ...,m k ) 



mil ■ ■ ■ m k \ 

x {[a + k - 2q + l) mi ■■■(a + k- q) mq - q +i) 



x 



((a + k-2q- m q+1 ] 



m q+ i+q 



la 



+ k-2q- m k ) mk+2q -k+i) , (4.66) 
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with a defined in Eq. (I3.39p . and prove it by induction over k. 

This conjecture is verified in the k = 1 and k = 2 cases. Since it is also valid for the two 
previous examples (A; — q — or 1), it only remains to prove it whenever k — q > 2. In such 
cases, the last two seed functions (pk-i and (fi k , used in Ath-order SUSYQM, are necessarily 
of type II. Suppose that Eq. (I4.66P is verified for k! < k. We may therefore write 

(«+2)/m - ( i W A ( m i> • • • » m q )A(m q+ i, m fc _ 2 ) 
x ((a + k - 2q + l) mi ■■■(a + k- q) mq - q+ i) 

x ((a + k - 2q - m q+1 ) mq+1+q ■ ■ ■ (at + k - 2q - m k „ 2 ) mk2+2q - k+3 ) , (4.67) 

(«+l)/W> - f i ^' A( ^ mi ' 111? m g )A(m g+ i, • • • , m k-i) 
x ((a + A - 2g + l) mi • • • (a + A; - g) mg _,+i) 

x ((a + A; - 2q - m q+1 ) mq+1+q ■■■(a + k-2q- m k ^ 1 ) mk _ 1+2q ^ k+2 ) , (4.68) 



and 



with 



(p+VfriS - i 1 ^ A ( mi > • • • » m q) A ( m q+i> • • • » m k-2,m k ) 

mi! • • ■m k _ 2 \m k \ 

x ((a + A; - 2g + • • • (a + A - g) m9 - f;+ i) 

x ((a + A - 2g - m g+ i) mg+1+g • • • (a + k - 2q - m k _ 2 ) mk2+2q _ k+3 ) , 

x (a + k - 2q - m k ) rrik+2q _ k+2 , (4.69) 



°' = Et g +i m, + g(A; - g - 1), 

^ = ESn'm + ?(fc-«-2) I (4-70) 

^' = Sf= g +i mi + m k + q(k - q - 1). 
From Eq. (14.571) . we obtain the recurrence relation 



#(0) = mk ;™ 1 k - 1 ti +1 \0) 9 £ +1 \0)/ 9 ^ 2 \0), (4.71) 
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which combined with the previous expressions yields 



A(mi, ■ ■ ■ , m q )A(m q+1 , m k _ 1 )A(m q+1 , . . . , m fc _ 2 , m fc ) 
A(m g+ i, . . . ,m fc _ 2 ) 
x ((a + fc - 2g + l) mi ■■■(a + k- q) mq -q+i) 

x ((a + k - 2q - m q+1 ) mq+1+q ■■■(a + k — 2q — m k _ 2 ) rnk ^ 2+2q „ k+3 ) 

x (a + k - 2g - m fe _.i) mfc _ 1+2g „ fc+2 (a + fc - 2g - m k ) mk+2q ^ k+2 . (4.72) 



Here we note that 

fc-2 



a' + a' — a" = + m fc „! + m k + q{k — q) = a, (4.73) 

l=q+l 

(a + k - 2q - m k ) mk+2q -k+2 , , , v (Af7AS 

— = (a + k - 2g - m fc ) mfc+2 g_ fc+1 , (4.74) 



a + 1 



and 



, . A(m q+1 , m k -i)A(m q+ i, m k _ 2 , m k ) 

(m k - m h -i) — s r 

A(m q+1 , . . . ,m fe _ 2 ) 

(fc-2 fc-1 \ / fc-2 

n n ( m i - m ) ( n (^ - m ^ 
i=g+lj=i+l / \i=g+l 

fc-1 fc 

= [ [ ] [ {mj - mj) = A(m 9+ i, . . . , m k ). (4.75) 

The final result for g^\o) is therefore given by Eq. ( 14.661) . showing that the conjecture is 
still verified for k, which ends its proof. 

Let us observe that in terms of a', Equation (I4.66P can be rewritten as 

(«}/m / y A(m 1 ,...,m q )A(m q+1 ,...,m k ) 

9f * { } { ' m l \---m k \ 

x ((«' + l) mi •••(«' + q) mq - q+ i) 

x ((«' - m q+1 ) mg+1+q •••(«'- m k ) mk+2q _ k+1 ) . (4.76) 

Since the angular momentum V in the starting potential (x) is necessarily nonnegative, 
the inequality a' > —1/2 is fulfilled. Furthermore, we have assumed < m\ < ■ ■ ■ < m q 
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and < m q+ i < ■ ■ ■ < rrik < a'. As a consequence, the two Vandermonde determinants in 
Eq. ( I4.76P are positive and 

(a + l) mi •••(«' + q)m q - q +i(a' - m q+1 ) mq+1+q •••(«'- m k ) mk+2q -k+i > 0. (4.77) 

We conclude that the sign of gjj*\o) is given by (— 1)°", hence is identical with the one of Cji* 
(see Eq. (I3.39j) ) . The combination of this result with the disconjugacy theorem of Sec. II 
ensures the regularity of the extended potential V^(x). 

V CONCLUSION 

In the present work, we have established some further properties of the function g^(z) 
appearing in the denominator of rationally-extended isotonic oscillator (or radial oscillator) 
potentials constructed in fcth-order SUSYQM [31J or using the multistep Darboux-Backlund 
transformation method [32]. We have first confirmed that \z) cannot have any pole at 
the origin and is a polynomial with non-zero constant term, a property that was conjectured 
in Ref. [31J. We have then determined the value of this constant term by induction over 
k and shown that it has the same sign as the coefficient of the previously determined 
highest-degree term z^ [31 J. 

These results have been combined with the disconjugacy theorem applied to Schrodinger 
equation [32] to prove the absence of zeros of g^\z) on the positive half-line and, therefore, 
the regularity of the corresponding rationally-extended potential. This has provided a 
demonstration of the latter property, which is more explicit than that already given in 
Ref. [32]. 

Our work has illustrated the power of the disconjugacy properties of second-order dif- 
ferential equations of Schrodinger type to check the regularity of rational extensions of 
well-known quantum potentials. Other examples of application might be found in cases 
where no direct information on the denominator zeros is available. This will be the object 
of further investigations. 
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